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Q ■ We propose to combine and slightly extend two existing "Les Houches Ac- 

cords" to provide a simple generic interface between beyond-the-standard- 
model parton-level and event-level generators. All relevant information — par- 
ticle content, quantum numbers of new states, masses, cross sections, parton- 
' level events, etc — is collected in one single file, which adheres to the Les 

Qh| Houches Event File (LHEF) standard. 

in: 



Abstract 



1. INTRODUCTION 



The simulation of interactions at the LHC is characterized by the use of many different programs spe- 
cializing in different stages of the calculation, such as matrix-element-level event generation, decay of 
resonances, parton showering, hadronization, and underlying event simulation. The communication of 
simulation parameters between those stages can be complicated and program-specific. For supersym- 
^ ' metric models, this situation has been greatly improved by the introduction of the SUSY Les Houches 

Accord [1] (SLHA) and its upcoming extension [2|. For general models however, there is still no cor- 
^ ! responding standard. In this note, we suggest an addition to the SLHA to allow for the specification of 

\ the quantum numbers, masses, and decays of arbitrary new states, thus generalizing the accord beyond 

O ■ its original supersymmetry-specific scope. We also make a proposal for how to include these model pa- 

rameter files into Les Houches Accord event files [3J (LHEF) in a standardized way. This both reduces 
K / \ the number of files that need to be passed around and minimizes the possibility for error by keeping all 

■ relevant model information together with the actual events. 

2. DEFINITION OF THE INTERFACE 

The concrete proposal consists of the following three points: 

1. Introduce new SLHA-like blocks QNUMBERS (for "quantum numbers") with the format: 

BLOCK QNUMBERS 7654321 # balleron 

1 # 3 times electric charge 

2 1 # number of spin states (23+1) 

3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 

4 # Particle/Antiparticle distinction {O=own anti) 
where this example pertains to a fictitious neutral spin-0 color-singlet self-conjugate particle to 
which we assign "PDG" code 7654321 and the name "balleron". That is, the BLOCK declaration 
should define a PDG code and, optionally, a human readable name after the # character (if no name 
is given, the PDG code may be used). We advise to choose PDG numbers in excess of 3 million 
for new states, to minimize the possibihty of conflict with already agreed-upon numbers [4|. The 



entries so far defined are: 1: the electric cliarge times 3 (so tliat most particles will have integer 
values, but real numbers should also be accepted); 2: the particle's number of spin states: 25* + 1; 
3: the colour representation of the particle, e.g., 1 for a singlet, 3 (-3) for a triplet (antitriplet), 8 
for an octet, etc.; 4: particle/antiparticle distincition, should be (zero) if the particle is its own 
antiparticle, or 1 otherwise. 

2. Use the existing SLHA blocks MASS and DECAY HI to define particle masses and decay tables. 
If the model in question is a SUSY model, a full SLHA spectrum [IJ can also be included. We 
propose that the reader should "turn on" SUSY whenever the SLHA SUSY model definition block 
MODSEL is present. 

3. Include the information from points 1 and 2 enclosed within the subtags <slha> </slha> in 
the <header> part of Les Houches event files ||3l. 

3. IMPLEMENTATIONS 

For the purpose of this contribution, the above proposal was tested explicitly by interfacing Mad- 
Graph/MadEvent with Pythia. Below we summarize the main aspects of these implementations. 

3.1 MadGraph/MadEvent implementation 

Starting from version 4 |5|, the multi-purpose MadGraph/MadEvent parton-level event generator 
by default includes a detailed summary of all simulation parameters in the output LHEF [3 | parton-level 
event file. From version 4.1.47, this information is stored in the XML <header> section. For the 
interface considered here, the relevant part of this section is a copy of the so-called param_card . dat 
MG/ME input file. 

The MG/ME param_card.dat uses an extension of the SUSY Les Houches Accord HI 13 
for model parameters in all implemented models. In particular, it always includes the SMINPUTS, 
MASS, and DECAY blocks. This file is used by MadGraph/MadEvent as an input for cross section 
computations and event generation but is not modified by the program. The file is instead assumed to be 
created by an external "Model Calculator". Such calculators are currently available on the web for the 
SM, MSSM and 2HDM models. Starting from the parameters in the Lagrangian (primary parameters), 
they calculate all needed secondary parameters (such as masses, decay widths, and auxiliary parameters). 
Note that widths and branching ratios can also be evaluated in an intermediate step by MG/ME itself or 
by external tools like DECAY or BRIDGE (6]. 

In previous versions of MadGraph/MadEvent, the param_card.dat file did not contain 
information regarding the particle content of the physical model considered. This information is stored 
in the particle, dat file filled by model writers during the model creation. Starting from version 
4.1.43, the template for inclusion of user defined models (called USRMOD) in MadGraph/MadEvent 
automatically generates the QNUMBERS blocks described above from the information contained in the 
particle . dat file. These blocks are then included in the default param_card. dat for the new 
model (and from there are copied into the LHEF output), such that no extra intervention is required 
to pass them to parton shower programs after parton-level event production. The script only outputs 
information for particles which have PDG numbers not identified as standard SM or MSSM particles, 
since those are assumed to be defined in the parton shower generators. 

Note that in the current version, the spin, color and particle/antiparticle information is automati- 
cally extracted, but not the electric charge, which is set to zero by default. This is due to the fact that, 
in MadGraph/MadEvent, the electric charge does not appear in the list of particle properties and 
is only defined through the value of the coupling to the photon. This issue will be addressed in future 
versions of USRMOD, but can currently be circumvented by fixing the electric charge information by hand 
at the end of the model implementation process. 



3.2 Pythia implementation 

The following capabilities are implemented in PYTHIA 6.414 fT\ and subsequent versions. 

Already for some time it has been possible to use the QNUMBERS blocks described above to define 
new particles in PYTHIA via its SLHA interface [8|. What is new is that, when reading an LHEF event 
file, Pythia now automatically searches for QNUMBERS blocks in the header part of the LHEF file, 
updating its internal particle data tables accordingly. It then proceeds to search for MASS and DECAY 
tables, and finally looks for other SLHA blocks contained in the header. If the SUSY model definition 
block MODSEL is found, SUSY is automatically switched on and the remaining SLHA blocks are read, 
without the user having to intervene. The read-in of LHEF files containing general BSM states, masses, 
and decay tables, should therefore now be relatively " plug-and-play". 

A note on decay tables: only 2- and 3-body decays can currently be handled consistently. They 
are then generated with flat phase space, according to the branching ratios input via the DECAY tables. 
The colour flow algorithms have been substantially generalized, but if too many coloured particles are 
involved (e.g., an octet decaying to three octets) PYTHIA will still not be able to guess which colour flow 
to use, leading to errors. Please also read the warnings in the section on decay tables in the SLHA report 
ITIl concerning the dangers of double counting partial widths and obliterating resonance shapes. To get 
around the restriction to flat phase space, either 1) use Pythia's internal resonance decays whenever 
possible (e.g., do not read in decay tables for particles for which Pythia's internal treatment is not 
desired modified), 2) perform the decays externally, before the event is handed to PYTHIA, or 3) do a 
post facto re-weighting of the generated events, based on the kinematics of the particle decays stored in 
the event record. 

The interfaces can of course still also be used stand-alone, independently of LHEF. The user must 
then manually open a spectrum file containing QNUMBERS and MASS information and give PYTHIA the 
logical unit number in IMSS (21) . New states can then be read in via either of the calls 

CALL PYSLHA (0, KF, IFAIL) ! look for QNUMBERS for PDG = KF 

CALL PYSLHA (0, 0, IFAIL) ! read In all QNUMBERS 

and MASS information can be read by 

CALL PYSLHA (5, KF, IFAIL) ! look for MASS entry for PDG = KF 

CALL PYSLHA (5, 0, IFAIL) ! read in all MASS entries 

where IFAIL is a standard return code, which is zero if everything went fine. (For read-in of a complete 
SLHA SUSY spectrum file, these direct calls should not be used, instead set IMSS (1 ) =11 before the 
call to PYINIT.) For stand-alone decay table read-in, the unit number of the SLHA decay table file 
should be given in IMSS (22) , and the corresponding read-in calls are 

CALL PYSLHA (2, KF, IFAIL) ! look for DECAY table for PDG = KF 

CALL PYSLHA (2, 0, IFAIL) ! read in all DECAY tables 



CONCLUSIONS AND OUTLOOK 

We have proposed a simple file-based interface between parton- and event-level generators focusing on 
the particular problems encountered in the simulation of beyond-the-standard-model collider physics. To 
deal with general BSM models, we add a new block QNUMBERS to the SLHA structure, which defines 
the SM quantum numbers of new states for use in subsequent resonance decay, parton showering, and 
hadronization programs. We also integrate the SLHA file into the existing LHEF format to minimize 
the number of separate files needed. The proposal has been tested explicitly by implementations in the 
MadGraph/MadEvent and Pythia6 Monte Carlo event generators. 

In the near future, also the Herwig-i-h lH and PythiaS HQ generators wiU be extended to auto- 
matically read in SLHA spectra from LHEF headers. Likewise, forthcoming versions of the CalcHEP 
ifTTT and CompHEP 1, 12 J parton-level generators will include write-out of this information in their LHEF 



output, including also the QNUMBERS extension. 

In the longer term, with the XML format emerging as the de facto standard for file-based interfaces, 
we note that it could be worth investigating the merits of formulating an XML-SLHA scheme, that is, 
transforming the current ASCII SLHA format conventions into a native XML form that could be parsed 
with standard XML packages. A concrete first realization of such a strategy is HepML |[T3]| which aims 
to unify the description of generator information in the form of standard XML schemes, in which an 
XML-SLHA scheme would form a natural part. The first release of the public HepML library has been 
implemented into CompHEP version 4.5, including also HepML headers in the LHEF output. 
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